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Abstract 

We show that the parton distribution functions (PDF) described by the statistical 
model have very interesting physical properties which help to understand the struc¬ 
ture of partons. The role of the quark helicity components is emphasized as they 
represent the building blocks of the PDF. In the model the sign of the polarized 
quarks PDF comes out in a quite natural way once the thermodynamical potentials 
with a given helicity are known. Introducing the concept of entropy we study the 
states made of|2n-|-d >, \u + d+s > and |2n-|-(i >, for a hxed Q^, the variation with 
X shows that the hrst state has a dominant entropy due to the effect of u quark. We 
prove that the PDF parameters obtained from experiments give in fact an optimal 
solution of an entropy equation subject to constraints. The same optimal property 
is proven for the structure functions and g^, and hnally to the quarks themselves. 
We develop a new approach of the polarized gluon density based on a neural model 
which explains its property, in particular, a large positivity value and an agreement 
with the positvity constraint. An extension of this neural approach is applied to 
quarks giving a coherent description of the partons structure. 
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1 Introduction 


The role of the parton distribution functions in QCD theory is essential to describe 
both unpolarized and polarized reactions, so a tremendous effort has been under¬ 
taken to hnd the most accurate distributions. In the litterature we remark that 
unpolarized and polarized PDF are treated as separated entities n-i. however, a 
simultaneous treatment of unpolarized and polarized PDF should in principle gives 
a more constrained determination. Since many years we have adopted this point of 
view in the framework of a statistical model where the PDF are built in from their 
helicity components. mm 

In the absence of a theory for the parton distributions two approaches are 
currently proposed. In one approach the distributions are approximated by different 
polynomials which require numerous parameters, with such an hypercube a carefull 
numerical analysis of the errors is necessary in order to get the most precise values 
of the distributions, however no attention is paid to the physical structure of the 
partons nor on the meaning of the parameters values. In our approach the physical 
structure of the distributions is introduced through a statistical model, this informa¬ 
tion allows us to work with a reduce number of parameters (21) whose meaning can 
be justihed. In both case a good description of the experimental data is obtained 
so we have a possible choice between a numerical formulation of the distributions 
versus a physical one where more emphasize is put on the partons structure [12] . 

The application of a statistical model, for instance, to a proton at rest which 
contains three quarks seems not justihed due to the low number of elements. But 
when accelerated in a collider the energy increase has not only an effect on its mass 
but also to create a large number of g g pairs or a quark gluon plasma which in a 
p-p collison materialize mainly in primary unstable particles observed in a detector 
as large number of tracks. These occurence of numerous pairs provide a justihcation 
for a statistical treatment of the partons interaction process. 

Let us mention an other application of the statistical model to different elastic 
scattering reactions in terms of quarks PDF dehned in impact parameter space, we 
have shown, in particular, that the gluon contribution is essential to explain the dip 
structure of the pp elastic differential cross section [13]. 

The paper is organized a follow: in section 2, the role of the thermodynamical 
potentials is discussed, in section 3 we consider the entropy of quarks states and show 
that the parameters values obtained from a ht correspond to a maximum entropy 
of these states. The same property is derived for the structure functions and 
Qp, and also for the quarks. In section 4 an analysis of the polarized gluon leads to 
dehne a neural model for its structure and in section 5 we developp an extension of 
this model to quarks. 
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2 The role of the thermodynamical potentials 


In the statistical model the thermodynamical helicity dependent potentials play 
an essential role in the construction of the polarized quark distributions and so have a 
direct consequence on the behavior of the polarized structure functions. The helicity 
decomposition of the quarks PDF is given by a quasi Fermi-Dirac distribution which 
is defined at the input scale = IGeV^ by the expressions: 



exp [(a: - X^)/x] + 1 


+ 


AqX^^ 

exp(a;/x) + 1 ’ 


( 1 ) 


xa^ix) = ^ _-_+_ ^ _. (2) 

X^ exp[(a: + XT)/a;] + 1 exp(a;/a;) + 1 

The last term is a diffractive contribution whose effect is to enhance the values 
of the unpolarized quarks at low x. The polarized quarks are dehned by the dif¬ 
ference xAq{x) = xq^{x) — xq~{x) and the unpolarized one by the sum xq{x) = 
xq^{x) +xq~{x), the antiquarks q have a similar dehnition. In these expressions x is 
a universal temperature, the introduction of the potential X^ and in front of 

Eqs. ([T]l2]) comes from the relation between the Transverse Momentum Distribution 
and the PDF [T3]. In a recent £t of unpolarized and polarised data made in Ref. 
[T^ we have obtained for the u, d, s potentials, 


= 0.475 ± 0.0008, X“ = 0.307 ± 0.001, = 0.245 ± 0.001, 

X" = 0.309 ±0.001, X+= 0.011 ±0.0008, = 0.015 ± 0.001. (3) 

From these values we obtain the potentials hierarchy 

X* < A',- < A+ < X„- < X- < A'+. (4) 


which is responsible of the quarks order of magnitude. Notice that the u and the d 
potentials are relatively stable since the analysis made in 2002 [10], it means that 
their values are a real intrinsic property of quarks, and they represent the master 
parameters of the statistical model. In Fig. [1] a plot of the polarized light quarks 
at = 4GeV^ is shown with their corresponding maximum potentials, we observe 
a correlation between the potentials values and the maximum or the minimum of 
the PDF. From Eq. ([T]) the sign of the polarized PDF Aq = q^ — q~ is related to 
the value of the thermodynamical potential helicity, more precisely, on the relative 
values of the potentials X+ and X~, the equality being excluded because it leads 
to a vanishing polarized PDF at the input scale, so we are led with two possibilities 
X~ > X+ and X~ < X+. In the case of u quark X~ < X+ so Au is positive, while 
for the d quark Xj > Xj" so Ad is negative, and for the strange quark s we have 
X~ > X+ leading also to a negative As (see HERMES experiment [T6]). 

^The PDF are evolved with the HOPPET program [15]. See Ref. [12] for more details. A 
Fortran program to compute the PDF independently of HOPPET is available upon request. 
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Figure 1: (color online) The polarized PDF at = 4GeV^ as a function of x. 
Vertical lines represent the maximum value of the potential for each quark. Shaded 
area uncertainty bands. 


For the antiquarks, the chiral structure of QCD gives an important prop¬ 
erty which allow to relate quark and antiquark distributions. The potential of an 
antiquark of helicity -h is opposite to the potential of a quark gf of helicity h 

. (5) 

So in the expression given by Eq. ([2]) the thermodynamical potentials have been 
interchanged with respect to the helicity ± and their sign taken to be opposite. This 
change of sign is due to the fact that a.q q pair can be created by a gluon through the 
process g ^ q + q. Due to the interchange of the potentials the sign of Ah is positive 
while Ad, As keep their negative sign. The respective signs are conhrmed by the 
parity violating asymmetry measured by the STAR polarized experiment [I^ 
in the process —)■ -|- X. 

Taking into account the numerical value of the potentials how they influence 
the spin structure functions, we will give two exemples. The polarized structure 
function xg{ has a maximum around x = 0.4, see Fig. [2], now we know that Au 
which gives the major contribution has a thermodynamical potential Xf^ = 0.46, 
so we observe a correspondance between this potential and the maximum in xg\. 
An other exemple is given by xg”, the data show that g^ is mainly negative over a 
large x region this fact can be explained by the inequality of the thermodynamical 
potentials X^^ < X^^ which implies A((i+ — d”) < 0, we also see in the figure that 
around x = 0.45 it has a postive maximum which reflects the influence of the u~^ 
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contribution at large x compared to Ad which is depressed in this region. 




Figure 2: (color online) The polarized structure functions xg^ (top) and xg'^ (bot¬ 
tom) as a function of x for Q'^ = 4GeV^. Experiments [l8]-[28]. Shaded area 
uncertainty bands. 
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Figure 3: (color online) The unpolarized PDF at = lOGeV^ as a function of x. 
Shaded area uncertainty bands. 


3 The quarks entropy 

In the previous section we have focused on the polarized PDF, the unpolarized ones 
are obtained from the relation q = q~^ + q~, we show in Fig. Elthe unpolarized PDF 
at = lOGeV^. With the parameters defined above they give a good description 
of the unpolarized structure functions in deep inelastic scattering, the neutrino cross 
sections, the neutral and charged current cross sections, and the jets production up 
to LHG energy, see Ref. [12]. 

We will explore a new property of the unpolarized PDF by considering a 
physical quantity precisely the entropy. The entropy can be calculated according to 
the definition given in Ref. [29] Eq. (16) 

E(Q^,x) = - ^ [xqi(Q^,x)\a(xqi(Q^,x)) + (1 - xqi(Q^,x))\n(l - xqi(Q^,x))] , 

( 6 ) 

where the sum runs over the quark components. We hrst remark that the vanishing 
of xqi{Q‘^,x) for x = 1, implies that E{Q‘^,x) = 0 in this limit. We propose to 
compute the entropy for the states made with \2u + d>, \u + d + s > and |2-u + J >, 
at a hxed = lOGeV^ as a function of x. In Fig. HI we see that the hrst state is 
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largely dominant over the last ones which seems to reflect the importance of matter 
over anti-matter. 

The curves shown in Figure 0] have been calculated with the values of the 
potentials obtained from a fit of experimental data discussed above, then a question 
arises, what is the origin of these values, does exist a possibility to obtain them 
independently of experimental data? In the PDF formulas described above we have 
introduced the following parameters: a normalization A, a power b of the variable x, 
a temperature x and the potentials. For a matter of simplification in the calculation 
let us assume that the following parameters A, b, x are held fixed to their actual 
values and consider now the potentials as free parameters which will be determined 
from a calculation of the optimal value of the entropy (E]) for a given value of x and 
Q^. In complete generality all the the parameters of the model should have been 
considered as free, but due the complexity of the computation we restrict our search 
only to the six potentials defined as the master parameters of the model. 

For this purpose we consider E{Q‘^,x) given by Eq. (ED as an objective func¬ 
tion which depends on u, d, s quarks subjects to the constraints 


0 < Xljg < 1 , 


Uy{x)dx = 2, 


dy{x)dx = 1, 


Au{x) > 0, Ad{x) < 0, As{x) < 0, 


(xu(x) -|- xd(x) -h xs{x))dx < 1, 


[s{x) — s{x))dx = 0 


( 7 ) 


The goal is to solve the sytem of equations (ED-Q with respect to the ther- 
modynalical potentials (supposed to be unknown) associated with u, d and s. The 
optimization is performed with the NLOPT software [30], which involves the objec¬ 
tive function, the constraints and their gradients with respect to the parameters. In 
addition, to confirm the results a brute-force method is also applied, it consists to 
And the maximum of the entropy by varying the parameters in a range defined as 
±50% of the fitted parameters values. 

We consider for a fixed = lOGeV^ a set of 20 x values in the range 
10“^ < a: < 1, the solutions for the optimal entropy are shown in Fig. |3]as circles 
for the state \2u d >, and squares for the state \u d s >, one observes that 
their values are close to the solids curves. These results show that the parameters 
obtained by this method have the same values (with an error around 2%) as the 
original ones obtained from a fit, so the entropy obtained from experimental data 
satisfies an optimal principle. 

We can envisage also to compute the entropy of a polarized state \2Au-\-Ad >, 
in this case there is a the difficulty which comes from the fact that to polarize, for 
instance, a proton one needs to apply a strong magnetic field, so there is coupling 
between the proton and an external field which introduces a complicated situation 
for the computation of the entropy because one has to disantangle the contribution 
coming from the external field and the other from the state itself. Nevertheless, 
we show in Fig. 0] the resulting entropy with a dash-dotted curve, the values are 
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Figure 4: (color online) Entropy at = lOGeV^ as a function of x for the states 
\2u + d >, \u + d + s >, \2u + d > and \2Au + Ad >, calculated with the experimental 
parameters. The optimal solutions of the entropy correspond to circles for \2u + d >, 
and squares for \u + d + s >. 



Figure 5: (color online) Difference of the entropy at = 1 and lOGev^. 
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much smaller than in the case of an unpolarized proton where the situation is more 
clear because the proton can be considered as in a free state. As a conclusion the 
calculation of the entropy obtained in an independent way from experiment has 
the consequence that the quarks PDF obtained from a £t correspond to a maxi¬ 
mum entropy principle, so the structure functions must share in some way the same 
property. We know that the entropy is sometimes associated with the disorder of a 
system and increases with energy. In Fig. [5] the difference of the entropy 6E between 
= 1 and lOGeV^ is effectively growing for the states discussed above. The state 
which involves the strange quark has the largest effect with respect to the disorder. 

A comparison of the optimum calculated entropy with experiment is not easy, 
an other test can be made with the structure functions. Using the same method 
as above, we consider again the thermodynamical potentials as free parameters in a 
certain range of values and search a maximum for the strcture functions Ff and gi 
for a fixed x and Q^. We hnd that the maximum is obtained when the potentials 
values are those obtained in the £t, so the following relation is derived 

maximum = F 2 £t = Uf experimental, (8) 

with the same relation for g\. 

In the above results the quark distributions are the essential source of infor¬ 
mation to obtain an optimum a property which should be reflected in the quarks 
themselves. To prove this we consider that in the unpolarized up and down quarks 
the potentials are now free parameters and a search is made for a maximum value 
given a fixed = lOGeV^ and x = 0.15. In Fig. [6] we plot the u and d values as 
a function of and X~ limited to a certain domain which generates a surface. 
Now, if we look for a maximum by imposing the constraints Eqs . [7] we obtain 
as a solution only one couple X^,X~ with a u and d values which correspond to 
those obtained in the fit (red point in the hgure). The same result is also obtained 
for the polarized Xu and Xd. So the optimum obtained for the entropy and the 
structure functions find its origin on the quarks properties. Let us mention that this 
maximum values of the light quarks unpolarized distriutions is also found with the 
parametrization MSTW 2008 |3T] and GT14[^ [32] . 

From this result we infer that nature tends to produce observable quanti¬ 
ties with a maximum probability taking into account some physical constraints, it 
remains to explain the origin of this effect. 


4 A neural model of polarized gluon 

The polarized gluon distribution is today not well known and subject to a large 
debate concerning its expression and sign. Our purpose is to clarify the choice made 
in our original model at the input scale and to propose a new interpretation in the 


thank J. Pumplin for informations on CT14 PDFs. 




Figure 6: (color online) Quark distributions as a function of X for = 
lOGeV^ and x = 0.15, left up-quark, right down quark. Maximum value red point. 


context of a neural structure. In the statistical model it is natural to assume a quasi 
Bose-Einstein distribution for both G and AG, so we dehne for the gluon at the 
input scale Ql = IGev^ 


xG{x,Ql) 


AqX^^ 

exp(x/x) — 1 ’ 


(9) 


and for the polarized distribution 


xAG{x,Ql) = Agx^°P{x) ■ 


1 

exp[x/x — 1) ’ 


where for P{x) we made the choice [5B] 


P{x) 


1 

(1 CgX^g'^ ’ 


( 10 ) 


( 11 ) 


notice that the introduction of an analogous rational multiplicative factor is also 
used in Ref. [3l]. A £t of polarized DIS data gives the values [12] 


Ag = 36.778 ± 0.085, = 1-020 ± 0.0014, Ag = 26.887 ± 0.050, 

Ig = 0.163 ± 0.005, cg = 0.006 ± 0.0005, dc = -6.072 ± 0.350 , (12) 


and for the temperature x = 0.090±0.002. We obtain a x^/d.o./. = 319/269. With 
these parameters AG{x) is positive a property conhrmed by experimental data from 
Hermes [55], Gompass [56] and A-j// from the STAR Gollaboration at BNL-RHIG 
[57] . In our analysis we take advantage that both the unpolarized and polarized 
gluon are determined in the same ht. We remark that unpolarized quarks and gluon 
are related through the evolution equations, in the same way the polarized quarks 
and gluon are related by an other set of evolution equations, now by construction 
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Figure 7: (color online) xAG{x, Q^) at = 2GeV^, dotted curve values multiplied 
by 10, lielicity components solid + and dashed curves as a function of x. Shaded 
area uncertainty bands. 



Figure 8 : (color online) The gluon helicity distribution xAG{x,Q‘^) versus x, at 
= IGeV^ (dash-dotted curve) and = 4GeV^ with cg 7 ^ 0 (solid curve), and 
cg = 0 (dashed curve). 
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Figure 9: (color online) Solid curve: the predicted double-helicity asymmetry 
with cg 7 ^ 0 for jet production at BNL-RHIC in the near-forward rapidity region, 
versus pt with the data points from Ref. HZI. 

Dot Dashed curve: the asymmetry with cg = 0. 


unpolarized and polarized quarks are related because the building blocks are the 
helicity components, as a consequence all the partons are linked together. 

A plot in Fig. [7] of the polarized gluon (dotted curve) at = 2GeV^ shows 
a maximum in the region x = 0.3, now from the values of G{x) and AG{x) we can 
deduce the helicity components (solid and dashed curves) their main difference is 
located in the same x region. 

A priori, it was natural to use for the polarized gluon an analogous expression 
like Eq. ([9]) for the gluon which is obtained by setting cg = 0 in P{x). With this 
assumption a new £t gives for the parameters 

Ag = (2.46 ± 0.1)10-^ be = 0.294 (13) 

with a x^/d-o.f. = 325/269, at this level we can consider that the 2 solutions fll2p - 
firSjl are equivalent, however, in Fig. | 8 ]we observe a marked difference for xAG. The 
pic obtained at = 4GeV^ when cg = 0.006 (solid curve) becomes a flat maximum 
when Cg = 0 and is reduced by a factor 4 (dashed curve). In order to separate the 
2 solutions we refer to the measurement of the double-helicity asymmetry for 
5 < Pt < 30GeV, in the near-forward rapidity region measured recently by the 
STAR Gollaboration [T7] . 

In Fig. [9] we have plotted the solution with cg 7 ^ 0 (solid curve )0 and the 
solution Cg = 0 (dashed curve). Despite large experimental errors the agreement 

^Results are taken from Ref. [33]. 
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with data is clearly in favor of cg 7 ^ 0 because the other curve grows too slowly with 
Pt, so the solution with cg = 0 is not acceptable. 

Let us now examine the function P(x) defined by Eq. (ITT]) . In Fig [10] a plot 
versus x at the input scale shows that P{x) is increasing with x, its hrst derivative 
is maximum for x = 0.41 and the second derivative (curvature) vanishes at the same 
X value. The shape of curve and the above properties show a close similarity with a 
sigmoid or logistic function whose basic expression is 

• (14) 

1 _|_ g AS 

A sigmoid function is used as an activation function in several domains, in particular, 
in neural networks applied to structure functions Pd, also in the exploration of 
opacity in elastic hadron scattering [ 8 ]. From our previous remark will consider that 
the P(x) function can now be replaced by a sigmoid of the form 

= 1 + 

where the parameter he dehnes a translation of the curve in the interval [0 < x < 1 ], 
so we dehne a new xAG at the input scale 

xAG'(x,Q 2 ) = ^( 2 :)—--. (16) 

exp(x/x)— 1 

A £t of polarized data yields the values 

A'g = 18.987 ± 1.5, cg = 77.94 ± 3.3, he = 22.18 ± 2.4 , (17) 

he is the same as in ((9|), we obtain a x^/d-o.f. = 326/271 very close to the original 

solution. 

For illustration we show in Fig. [TTjthe function S'(x). calculated with the 
fitted parameters flT7|) . It is interesting to compare xAG obtained with P(x) and 
cg 7 ^ 0 with the case of S'(x). In Fig. [12] the solution with P(x) corresponds to 
the dashed curve and with S'(x) solid curve, we observe that the latter has a more 
pronounced peak which decreases with and moves slowly toward smaller x values. 

The validity of the new polarized gluon can be tested by computing the 
asymmetry A-}//, the Fig. [13] shows a good agreement between the two solutions, 
P(x), Cg 7 ^ 0 and S'(x). 

In this section we have explored three possibilities to describe the helicity of 
the gluon. Starting with the expression Eq. f[TT]) used in [T2] which was phenomel- 
ogical, now we have shown that a more physical expression given by a sigmoid Eq. 
f[r5|) gives also a good desciption of polarized experimental data. 

How we can interpret the role of the sigmoid S (x). When 2 protons collide the 
gluon receives different fractions of the momentum coming from the quarks which 
are collected statistically with a Bose-Einstein distribution, next the function S{x) 
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Figure 10: (color online) The function P{x) lEq. fllip ) versus x for the polarized 
gluon. 



Figure 11: (color online) The activation function S{x) (Eq. (IT^ l versus x for the 
polarized gluon . 
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Figure 12: (color online) The gluon helicity distribution xAG{x, Q^) versus x, at 
= 4GeV^ calculated with P{x) (dashed curve) and with S{x) (solid curve). 



Figure 13: (color online) Dot dashed curve: Our predicted double-helicity asymmetry 
^LL with P{x), cg 7 ^ 0 for jet production at BNL-RHIC in the near-forward rapidity 
region, versus pt and the data points from Ref. ra- 
Solid curve: the asymmetry calculated with the function S{x). 


O Output A G(x) 

_ 

Figure 14: (color online) A schematic representation of a polarized gluon in a neural 
model. B — E{x) is the Bose-Einstein distribution, S{x) the activation function. 
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plays the role of an activation fnnction which synthesizes in an ontpnt signal AG. 
This mechanism allow us to define a neural representation of the polarized gluon 
whose schematic view is given in Fig. [TTl From this result it is tempting to apply 
the same representation to the unpolarized gluon, we now introduce in the gluon 
distribution an activation function fllSp . A global fit gives the parameters 

Ag = 33.23 ± 1.2, ec = 281.16 ± 6.1, he = -1.923 ± 0.01. (18) 

The resulting S(x) function is plotted in Fig. [T51 we observe for x > 10“^ a value of 
S(x) around 0.9-1 almost independent of x. It implies that the activation function 
makes no modiheation on the output distribution G, which can be interpreted by 
the fact that in order to maintain the conhnement of quarks any momentum transfer 
is allowed, also as stated in the introduction the creation of a maximum of qq pairs 
with increasing energy implies no selection. 



X 


Figure 15: (color online) The activation function S(x) versus x for the unpolarized 
gluon. 


In this new approach of the polarized gluon we would like to examine the 
ratio AG{x, Q^)/G{x, Q^) discussed in Ref. |33]. In Fig. [T6]the ratio is plotted as a 
function of x for four values, in this range the positivity condition AG{x)/G{x) < 
1 is statisfied, and for a hxed it increases with x, near the limit x = 1 the values 
are very close to 0.5. At the input scale the Bose-Einstein function cancels in Eqs. flU]) - 
flTB]) it results that AG/G{x = 1) —)■ Ag5G(l)/AGS'AG(l); now taking into account 
the ratio of the normalization factors we obtain the value 0.57, so the limit of the 
ratio AG/G is different from 1 as required by the counting sum rule. 


15 



1 


cr 

'o 10 “^ 


O' 

o 

<1 


10 


10' 




>-3 



10 


10 ' 


>-5 


10 ' 


>-3 


X 


Figure 16: (color online) AG{x,Q'^)/G{x,Q‘^) at = 2,4,10, lOOGeV^ (solid, dot¬ 
ted, dash-dotted, and dashed curves respectively) as a function of x. 

5 A neural model applied to quarks 

In the previous section we have focused on the structure of the gluon in a neural 
model, now a question arises for the quarks, can they share the same structure? 
The unpolarized quarks PDF are known with a good precision, and most of the 
parametrizations agree to produce the same values in x and Q^, in the polarized case 
there are more uncertainties but the observed shapes are more or less indentical, so 
the neural description we give in the gluon case seems not necessary. Nevertheless, 
looking at our PDF expressions Eqs. ([H [2]) we have the product of a Fermi-Dirac 
distribution by an helicity dependent funcion Ax^X^ for quarks and Ax^/X^ for 
antiquarks, so we can try to apply the same approach where the incoming momen¬ 
tum is collected now by mean of a Ferm-Dirac distribution and then filtered by an 
activation function to produce the quark distribution. Our objective is to obtain a 
coherent neural structure for all the unpolarized and polarized PDF. Several possi¬ 
bilities exist to introduce an activation function, we made the following choice where 
the original parton expressions for q, q and G are preserved when S{x) = 1. 


xq'^{x) = Sq{x) 



(19) 


xq^{x) = Sg{ 



—- ■ -h 

X^ exp[(a;-F XT)/x]-F 1 


exp{x/x) -|- 1 



( 20 ) 



( 21 ) 
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( 22 ) 


xAG'(x,Qo) = Sag{x)Agx’'° 
with an activation function defined by: 


exp{x/x — 1 ) ’ 


Si{x) 


1 

_|_ Q-eiX+hi ’ 


(23) 


the index i = q,q, G, AG, where Sg is identical for u, d quarks, and Sq for u, d To 
determine the parameters a global £t at NLO with an activation function included 
in all PDF is made with the same data set used in [12], we obtain a x^/d-o.f. = 
2506/2128 = 1.18. The thermodynamical potentials are slightly modihed 


A+= 0.540 ± 0.0014, X-= 0.336 ±0.0012, X+= 0.268 ± 0.0013, 

X" = 0.349 ± 0.001, X+= 0.0111 ±0.0011, X 7 = 0.0147 ± 0.0012 .(24) 


The activation function parameters are given in Table [T] and the corresponding 
functions Sg are shown in Fig. [T71 The curves characterize the response of partons 


i 


hi 

u, d 

27.16 ± 1.3 

0.7 (fixed) 

it, d 

23.37 ± 1.07 

D 

s 

15.27 ± 0.9 

D 

s 

8.34 ± 0.5 

D 

G 

281.67 ± 3.9 

-1.82 ± 0.1 

AG 

77.71 ± 2.0 

22.07 ± 1.24 


Table 1: Parton parameters of the activation functions. 


to a signal, the momentum, we observe a hierarchy where the u, d quarks have 
the dominant effect followed by antiquarks, strange and antistrange, it corresponds 
to the observed relative size of the PDF. In this first approach we have used the 
same activation function for u, d, Au, Ad, idem for the antiquarks, the strange 
and antistrange, but a more rehned version could introduce an activation function 
for each helicity components. 
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Figure 17: (color online) The activation function S{x) versus x for g, g, s, s, in the 



Figure 18: (color online) Unpolarized u, u solid curves, when 5 = 1 dashed curves. 
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Figure 19: (color online) Polarized Au, AG solid curves, when 5 = 1 dashed curves. 


In order to show more precisely the effect of the activation function we plot 
in Figs. [TSlIT^ at = lOGeV^ the PDF when Sq ^ 1, or when we set abruptly 
Sq = the effect of the activation function reduces their values at small x because 
at large x it becomes close to 1, the most striking effect is observed for the polarized 
gluon. 

We conclude that the model of a neural structure for the PDF is perfectly 
compatible with unpolarized and polarized experimental data. 

6 Conclusion 

The statistical model provides a better knowledge of the nature of the parton dis¬ 
tribution functions in the sense that their usual properties appears as a simple 
consequence of the statistical functions (Fermi or Bose-Eistein) and the thermody¬ 
namical potentials. The model gives a fairly good description of unpolarized and 
polarized experimental data with a reduced number of parameters and also presents 
a good laboratory to explore the partons structure. The sign of the polarized PDF 
for the quarks is hxed by the potentials and the dominance of the unpolarized and 
polarized u over the d appears in a quite natural way. The calculation of the entropy 
for the two states \2u + d > and \u + d + s > satisfies a maximum entropy principle 
with the potentials obtained from the experimental value of the PDF parameters. 
We have also proven that this optimum principle is valid for the structure functions 
Fp, Qp and at the end to the quarks themselves. 

A description of the polarized gluon in term of a neural model gives a more 
physical insight on its strucrure and removes the arbitrariness of the original for- 
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mulation. An extension of the nenral approach to qnarks is derived leading to a 
coherent pictnre of the partons structure which describes both unpolarized and po¬ 
larized experimental data. From a pure numerical point of view the polynomial and 
the statistical approaches give the same results, however, the last one provides a 
new explanation of the parton structure. It is clear that a neuron is not a parton 
but we have shown that the mathematical formulation applied to the former can be 
extended to the later. This hrst approach certainly needs further developments by 
considering helicity dependent activation functions, and also an extension to heavy 
quarks has to be envisaged. 
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